Quasiclassical Born-Oppenheimer approximations 
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We discuss several problems in quasiclassical physics for which approximate solutions were recently 
obtained by a new method, and which can also be solved by novel versions of the Born-Oppenheimer 
approximation. These cases include the so-called bouncing ball modes, low angular momentum states 
in perturbed circular billiards, resonant states in perturbed rectangular billiards, and whispering 
gallery modes. Some rare, special eigenstates, concentrated close to the edge or along a diagonal of 
a nearly rectangular billiard are found. This kind of state has apparently previously escaped notice. 



I. INTRODUCTION 

A major success of the quasiclassical method began 
with Martin Gutzwiller. With his famous 'trace formula' 
he started on the road to quantize 'hard' chaotic sys- 
tems, ones possessing only strongly unstable orbits. The 
work he did and the work he inspired can now be said 
to have solved the quantum spectral problem for two di- 
mensional hard chaos systems. 

The study of Gutzwiller's problem led to the introduc- 
tion of new tools such as the dynamical zeta function J2| . 
The surface of section transfer operator (SSTO), a gen- 
eralization of the boundary integral method, is an impor- 
tant technique invented by Bogomolny S. This operator, 
T(E), is a Fredholm kernel which carries a ray of energy 
E quasiclassically from one intersection with a surface of 
section to the next. It was first employed as a power- 
ful reformulation of the Gutzwiller trace formula. The 
spectrum is given by the zeroes of the Fredholm deter- 
minant D(E) — det [1 — T(E)\ which is an expression of 
the dynamical zeta function pi . The trace formula itself 
is given by the logarithmic derivative of D(E), expanded 
in traces of powers of T. These traces are expressed qua- 
siclassically in terms of periodic orbits. 

However, chaos is not required to use the SSTO. It 
was discovered |jj, || that sometimes the equation 



T(E)ip = vb 



(1) 



could be solved quasiclassically for both the energy 
[equivalent to D(E) = 0] and the 'surface of section' 
[SS] wavefunction vb. The two-dimensional wavefunction 
is obtained from vb. In this way we found wavefunctions 
and energy levels for some problems that had been exten- 
sively studied by numerics and trace formulas M, M, || , 
]9| , but which had not been suspected of having simple 
analytic solutions. After obtaining these solutions, we 
realized that they could often also be found to the same 
level of accuracy by extensions of the Born-Oppenheimer 
approximation (BOA) poft . The BOA is not usually re- 
garded as quasiclassical, of course. Further research re- 
vealed that these extensions of the BOA are related to 
the 'parabolic equation' and etalon methods mentioned 
below. This paper is devoted to illustrating these rela- 
tions. 



The problems susceptible to this method are ones 
which are 'locally' 'nearly integrable'. By this we mean 
that there is, for these problems, a region of phase space 
where the short orbits are close to those of some inte- 
grable system. This often gives rise to rare special states, 
with striking properties. Perhaps the best known exam- 
ple is the 'whispering gallery' idea of Lord Rayleigh |ll| ] 
which explains quantum phenomena 113] as well as effects 
in acoustics and seismics, optics ]13[ and radio propaga- 
tion. 

A number of methods have been developed to deal 
with such states, where it is as important to understand 
the wavefunctions as it is the energies. Among these 
methods are the Keller- Rubinow or ray method p3, the 
parabolic equation method (PEM) of Leontovich Ma and 
Fock |l6[ , and the etalon method of Babich and Buldyrev 
|p7| . However, there is no method so systematic that it is 
sure that all such special states have already been found. 

The Keller-Rubinow method is purely quasiclassical 
and is based on first understanding the classical mechan- 
ics, or 'rays' in the optics nomenclature. It sorts out 
the local integrability, exploiting caustics or an adiabatic 
invariant of the motion. The last two methods men- 
tioned start from the basic partial differential equation 
[e.g. Schrodinger or Helmholtz] defining the problem, 
which is rescaled in appropriate variables and approxi- 
mated. Systematic corrections to the leading order are 
also studied, and indeed seem to be the focus of much 
of the theory. These latter methods are related to the 
BOA, a fact not previously noted. 



II. BORN-OPPENHEIMER APPROXIMATION 



The Born-Oppenheimer or adiabatic approximation is 
fundamental to the quantum theory of molecules as well 
as to the quantum theory of the solid state. It treats the 
electrons' position r e as 'fast' variables , compared to 
the 'slow' ionic positions R^. This is based on the small 
electron-ion mass ratio m e /Mi. 

Most texts give a simple formulation of the BOA p8| . 
Schrodinger's equation is, in terms of electronic and ionic 
positions r e , Rj, 



n z 



V 2 + V(r e ,R,) 



2m e c v ' " 2M, 
= ^(r e ,R l ). 

The Born-Oppcnhcimcr Ansatz is 

*(r e ,R 4 ) = $(r e |R l )V'(R l 



*(r e ,Ri) 



(2) 



(3) 



where $ is, say, the N'th eigenstate in the electronic 
variables which solves 



2m. 



-V= + V(r e ,Ri) 



$(r e |Ri) = U(Ri)$(i e \Ri). 



(4) 

In $, Ri is treated as a parameter, i.e., the quantum mo- 
tion of the fast variable r e is found for a fixed value of the 
slow variable Ri. The 'potential' U is an energy eigen- 
value of this 'fast' part of the Schrodinger equation, often 
but not necessarily that of the electron ground state. The 
adiabatic invariance is invoked by the assumption that 
the electrons' quantum state label N does not change as 
Ri is slowly varied. 

Now use Eq. (g) in Eq. (g), neglect derivatives of <I> 
with respect to Rj, multiply by $(r e |Rj), and integrate 
over r e . This gives for the slow variables Ri the equation 



2Mi 



U(K t ) 



<A(Ri) = £ty(R0, 



(5) 



which completes the leading order solution. 



III. IMPERFECT SQUARE CAVITY 

As a first almost trivial example with apparently new 
results, consider a two-dimensional microwave or laser or 
quantum-dot cavity which is not quite square. For ex- 
ample, with Dirichlct's conditions, take a trapezoid with 
boundaries x — 0, x = L, y = L, y = ex where e is 
small. Consider states whose wavelength A is short, i.e. 
A << L. Even if the perturbation is so small as to sat- 
isfy eL << A, there can be states which differ drastically 
from the states of a perfect square. The condition for 
such novelty is y/eL > A. 

The trapezoid problem was studied formally in Morse 
and Feshbach fl9| , and is also the basis of recent 
work on quantum chaos ||, but no special states are 
mentioned. Look for a state fast in the y direction, 
slow in x. Take L = 1, and choose for the Ansatz, 
$(t/|x) ex sin [nn (1 — y) / (1 — ex)] with mr large, satis- 
fying d 2 <fr/dy 2 = 'constant' x$ and vanishing at y = 1 
and y = ex. The 'constant' is of course a relatively slowly 
varying function of x. Then U(x) — [rnr/ (1 — ex)] ~ 
?i 2 7r 2 + 2en 2 TT 2 x. The resulting equation for ip is familiar, 



ip" + ofxip = E m ip 



(6) 



where E m = E — n 2 i: 2 and a — (2en 2 -K 2 ) 1 ^ 3 . The solution 
is (if a >> 1) an Airy function, ip(x) = Ai(air — z m ), 



where E m — a 2 z m . Let Ai(— z m ) = 0, to make ^(0) = 0. 
If a — z m > > 1 , ip also effectively vanishes at x — 1 . 

Thus we have a number of rather simple special wave- 
functions, labelled by n, to, which are concentrated along 
the long edge of the trapezoid. Namely, 



*, 



sin [nir (1 — y) / (1 — ex)] Ai(aa; — z m ) (7) 



We show in Figs. 1,2 a couple of representations of these 
states, for n = 55, m = 1, 2. We choose e = 0.01, rather 
less than A — .036, which in turn is rather less than 
^fe = 0.1. The wave functions were obtained numerically 
and are compared with Eq. (m). Presumably a result like 
this, so fundamental to the construction of cavities, is in 
an engineering textbook somewhere, but if so, we haven't 
found it. 




x,y 

FIG. 1. Numerically obtained states for the nearly square 
trapezoidal billiard, e = .01, n — 55, m = 1,2. a) 
*55,i(a; = 0.01,?/) vs y; b) * 55 ,i(a;,y = .99) and Ai(ax - z±) 
vs x; c) *55,2(£, y — .99) and Ai(ax — Z2) vs x. The magnitude 
of the states has been normalized. 

It is not necessary to have an explicit small parameter. 
Keller |14| obtained approximate energies for 'bouncing 
ball' states in convex billiards. Even more dramatic re- 
sults were obtained by the BOA some time ago [^(J ex- 
plaining the 'bouncing ball' modes observed numerically 
|2lf in the Bunimovich stadium billiard. These results 
may be readily be generalized to a distorted stadium. 
Following Primack and Smilansky Q, take a 'slanted' 
stadium billiard in which the radius of the 'left' semicir- 
cle is R{1 - e/2), and the right radius is R(l + e/2). The 
two endcaps are separated by a distance 2a. The BOA 



.y+R-jjx) 



Here y = R — £,(x) is the 
This gives a potential 
7r 2 /2i? 3 ) £0). For the 



is $(y\x) = sin r » 2(H _ c(g)) 

upper boundary of the billiard 

for the slow equation V(x) = ( 

standard stadium, £ = for \x\ < a and becomes positive 

outside that region. If n 2 ir 2 is large, V(x) rises rapidly 

for |a;| > a, so it is a sort of square well potential. For the 

slanted stadium, the bottom of the well V(x) is sloped, 



m 



cR/a fo 



r a; < a. 



and again the poten- 



tial rises rapidly outside that region. The effect of this 



slope depends on a/R and the slow quantum number m 
as well as n. Thus, if the change of 'potential energy' 
V(a) — V(— a) ex (n 2 / R 2 ) e >> (m/a 2 ) , the spacing of 
the slow energy levels, there will be a large effect on the 
slow wavefunctions, concentrating them in the wider part 
of the billiard, similar to the trapezoidal billiard above. 
With the opposite inequality, the slanting of the sides of 
the billiard can be neglected. 



the phases which appear in the transfer operator and its 
eigenf unctions according to their rate of change as a func- 
tion of position on the surface of section. Thus, there are 
'fast' and 'slow' parts to the Bogomolny equation. 

In the trapezoid problem above, we take the space 
part of the surface of section to be y — 1. The SSTO for 
billiards is 



T(x,x'\k) 



1 


d 2 S(x,x') 


2ixih 


dxdx' 



.■S(*,*') 



(8) 




FIG. 2. Density plot of ^55,2 {x, y) | , the same state as in 
Fig. lc. The state is shown in the lower quarter of the billiard, 
y < 0.25, and the scale is expanded in the y direction by a 
factor 3 in order to show more details of the wave function. 
A dashed line is shown at y = 0. 

This smooth transition of the quantum levels contrasts 
with the mathematics of the classical periodic orbits used 
in the trace formula. With parallel sides to the billiard, 
there is a set of nonisolated periodic orbits contributing a 
very large term to the trace formula. An arbitrarily small 
slope mathematically eliminates all such nonisolated pe- 
riodic orbits. However, for a bounce or two, they remain 
close to the ideal periodic orbits of the channel with two 
parallel sides. 

Reference M sorts this out with respect to the trace 
formula. It is remarkable that the eigenstates and 
energies start to be substantially modified when e > 
(R/amr) while the low terms of the trace formula and 
thus the so-called length spectrum, are only modified if 
e > R/anir. In other words, the low terms of the trace 
formula are essentially unmodified if the actual billiard is 
a ideal stadium billiard with 'optically flat' errors much 
smaller than a wavelength. This is however, not good 
enough to guarantee that the wavefunctions are close to 
the ideal wavefunctions. 



IV. BOGOMOLNY OPERATOR 

The preceding results can be obtained by a method 0] 
based on Bogomolny's |3| quasiclassical surface of sec- 
tion transfer operator, T. The main idea is to organize 



Here the action, S(x,x') = hkL(x,x'), is that of the 
classical orbit leaving the surface of section at x' and 
returning to it for the first time at x. It is expressed by 
L(x,x'), the length of this path. The wavenumber k = 
\JE in our units. We focus on orbits in the neighbor- 
hood of the (1,0) resonance orbits of the unperturbed 
square. Such an orbit starts at r'= (x', 1) and returns 
to r = {x, 1) after making one bounce from the bottom 
at approximately r = (5(0; + x'), \e{x + x')). For small 
x - x' , L ?x 2 + j(x - x') 2 - e(x + x'). 

The second term of kL is rapidly varying with x' , 
the last term is relatively slow. We solve Eq. (|l]), i. 
c. J dx'T(x,x'\k)i/j(x') — ip(x), by making the Ansatz 
ip{x') = exp[±ik-\/eF(x')}. This function, which has 
an intermediate rate of variation, followed by station- 
ary phase approximation for the integral, causes the 
important values of x' to be close to x, and allows 
F(x') rs F(x) + {x 1 — x)f{x). Solution of Eq. (|l|) requires 

ky/eF — 3 [z m — ax] , where a is defined in the previ- 
ous section. This gives for ip the WKB approximation to 
the Airy function. 

In short, the simplest version of this SSTO tech- 
nique gives quasiclassically the same ip as the BOA. 
Knowing ip, the full wave function is given by 
^(r) = J dxG [r,r(x)]tp(x), where the kernel G is re- 
lated to the free space Green's function between a point 
r in the interior of the billiard and the point r(x) on the 
surface of section. The result is essentially the BOA. 

Note that the SSTO gives first the surface of sec- 
tion wavefunction, then the BOA wavefunction, reversing 
the order of the Born-Oppenheimer approximation. The 
BOA is simpler and more intuitive, when it works. 



V. DIAGONAL STATES 

We next give a result for which we cannot construct 
a simple BOA Ansatz. Consider the same trapezoidal 
billiard but now look for special states related to the (1,1) 
periodic orbits of the square. Then, it's not too hard 
to work out the SSTO by making the Ansatz ip(x') — 
exp(±ikx' J \/2) exp [±iky / eF(x')]. The first phase factor 
chooses orbits at approximately 45° to the square sides. 
In this case, the result is that there are states, induced 
by the small perturbation, that are concentrated along 
the diagonal (0, 0) -> (1, 1) but not along (0, 1) -> (1, 0). 
[In an appropriate parameter range, we do expect states 
with 'scars' along that diagonal, however.] Such a state, 



obtained numerically, is shown in Figs. 3 and 4. We 
have found the 2-dimensional state ^(x,y) to order e, 
but it is too cumbersome to display here. The following 
wave function (for even n) vanishes on the boundary, and 
is the same as our theoretical wave function to order e° 
(but not e). It is not an adequate BOA Ansatz but it 
is simple and it gives the gross features of the result. 
The function is ^>(x,y) *=s cos [iri7r(a; + y)~\ u(x — y) — 
cos [^n7r(a; — y)] u(x + y) where u[x) = Ai(a \x\ — Zq) for 
\x\ < 1, extended with period two outside this range, 
and y = (y — ex)/(l — ex). Here zq ~ 1-02 is the first 
maximum of the Airy function. The energies of such 
states are given by E n _ m = i (wr) + ^a 2 z m where z m is 
a zero or extremum of the Airy function. We display the 
remarkable nodal structure of such a state in Fig. 5. 



|3<x,y) 




FIG. 3. A three-dimensional representation of a state of 
the trapezoid concentrated near the long diagonal. The the- 
oretical state is plotted, e = 0.03, n = 50, m — 0. 



angle. First, we consider an auxiliary problem, a method 
of images, from whose solutions the desired answer is 
constructed. 

The auxiliary problem extends the domain of the wave- 
function to all of two dimensional space. The perturbing 
magnetic flux, originally defined only in the rectangle, 
is reflected about each of the sides of the rectangle, and 
then the process is repeated so that the result forms a 
lattice of period 2a in the x-direction, and period 2b in 
the y-direction. The vector potential is chosen to have 
this periodic structure. 




With a little practice, it is easy to foresee when special 
states exist or not. For example, there are no states con- 
centrated along the diagonal for the symmetric trapezoid 
billiard with sides x = 0, 1, y = ex, 1 — ex, but there are 
states concentrated along the side x — 0. For the paral- 
lelogram billiard with sides x — 0, 1, y — ex, y = 1 + ex, 
there are states along the long diagonal but no states 
near the edges. 



VI. EXTENDED ZONE SCHEME 

The BOA requires a choice of fast and slow variables. 
In the preceding example, a reflection of the wave inci- 
dent on a boundary at 45° switches the roles of x and y. 
The following device overcomes this difficulty, although it 
is cumbersome to apply to the example of the preceding 
section. 

Consider a perturbed rectangular billiard. The per- 
turbation can be a potential, a magnetic field, or mag- 
netic flux lines. We assume the perturbation is classically 
weak, meaning that the classical orbit does not change 
much because of the perturbation in one traversal of the 
system. Our SSTO technique automatically gives the re- 
sult |2^], where pictures of some special states may be 
found. 

For the BOA, we choose the important example where 
a magnetic field perturbs the (1, 1) resonance of the rect- 



FIG. 4. The same state as shown in Fig. 3. The numeri- 
cally obtained state is plotted, basically along the diagonals 
of the trapezoid, a) ty(x,x), b) ^(x, 1 — x). Also shown in 
part b) is the theoretical wavefunction, which is smoother. 
Part b) is plotted for < x < 1/(1 + e). 

Next we find the high energy wavefunctions of a peri- 
odic lattice, where the lattice potential is classically weak. 
The methods developed in the study of channelling |23] 
are appropriate and are effectively the BOA. We start 
with the (1, 1) channel, which amounts to assuming that 
the fast direction is given by the variable £ = (ax + by)/c 
and the slow by r\ — (bx — ay)/ c, where c = y 'a 2 + b 2 . 
Take charge e = 1, mass = i, and h — 1 and rewrite 
Schrodinger's equation in variables £, r\ as 

{[-»$ - A^i, V )f + [-id n A n (t, v)} 2 } *& V) 
= E*(Z,V)- (9) 

The channelling approximation averages along the fast, £, 
direction. This amounts to assuming $(£|t?) = exp(i/c£) 

and replacing A 5 (£,7y) by I e (r/) = /„ e d£A e (£,r?)/L s . 
Here L^ — 2c. Because k = —id^ is large, the terms 
can be neglected in comparison to kA^(ri) and 
The remaining rj dependence is obviously 
periodic. 
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FIG. 5. A density plot of a state similar to that shown in 
Figs. 3,4, with n = 28. The dashed line is y — 0. 

Take the origin at the rectangle center, for symmetry 
reasons. Let ip{r)) = u{cq/b — a/2) = u(x — ay/b — a/2). 
Then u satisfies 



u'^{x) + V{x)u m {x) = E m u m {x) 



where 



V(x) 



_2*> l> A ( bx 



rkAt (% 



2c) 



(10) 



(11) 



The reason for the shift in the definition of u is that then 
V{x) — V(—x), and u will be simple on the boundaries 
where the Dirichlet conditions are imposed. 

The integral denning A^ is the flux enclosed by the pe- 
riodic orbit in the original rectangle, and so is indepen- 
dent of gauge. For uniform field B, A^(r/) = +Brj(l — 
cq/ab) for < r\ < ab/c, A^(i]) = —B(2ab/c — r))(cr)/ab — 
1) for ab/c < r/ < 2ab/c, A^(r] + 2ab/c) — A^(r/). In terms 
of x, V(x) = -{Bkb 3 a/c 3 ) [i - 2(a;/a) 2 ] for |ar| < \a, 
repeated antiperiodically outside this region. The solu- 
tion of Eq. ( |l0|) is a Bloch state, which satisfies 



u m {x + 2a) = e 1 ' u m (x) 



(12) 



where to labels the 'band' and /3, the 'crystal momen- 
tum'. The band index to is a quantum number and f3 is 
yet to be determined. The energy is given by 



E 



c 2 E„ 



jb\ 



(13) 



Quantum states will be strongly affected by the flux 
if the dimensionless parameter Bk(ab/c) 3 is large. This 
parameter is 2n((j)/(f>o)(kc) [ab/c 2 ) , where <j) is the flux 
4> = abB, and 0o = 2i:h/e is the flux quantum. 

We therefore find a set of approximate solutions to 
the plane problem, labeled by k, to and /3, ^j(x,y) — 
exp [ik{ax + by)/c] u m (x—ay/b—a/2). Given \P/, another 
solution is given by the symmetry under rotation by 7r, 
$ iu(x,y) = *&i(— x, — y). Two other solutions are given 
by ^n(x,y) = exp [ik(— ax + by)/c] u m (—x — ay/b + a/2) 
and *i>iv(x,y) = ^n(—x, —y). All four of these channels 
clearly have the same energy. 

The solution for the original rectangle can be con- 
structed from these four solutions, provided some quan- 
tization conditions are met. We consider two subspaces, 



r = ±1, even or odd under (x, y) 
can be written in the form 



-y). The state 



*, 



A (*/ + r*///) + B(^ n + r^ IV ). 



(14) 



Imposing the condition \l/ r , 
- B exp(-ixa) , and exp(i/3) 



[±\a,y] = yields A = 
exp(-2ixa) where Xa 



ka 2 /c. The condition ^> nm {x,±\b) = gives A = 
—rBexp(ixb) an d exp(i/3) = exp(2ixb) where Xb — 
kb 2 /c. Therefore, Xa + Xb — ^n, where n is even for 
r = 1, and odd for r = — 1. This yields 



Tin 



(15) 



One finds f3 = 2xb — 2-kth, — ~2x, a + 2irn a . We may 
assume that \/3\ < n, which fixes the integers n a .b as the 
integer part of (xa.b^ 1 + 1/2). Then n a + n^ = ck/n = 
n. 

For a square, a = b, there is a further symmetry, and j3 
is either or w, that is, the appropriate state is either at 
the band top or band bottom. For an arbitrary ratio a/b, 
k£ = k(ax + by)/c is not of the form Ttn a x/a + Tin^y jb. 
However, the n a .b just found are the integers which most 
closely satisfy this relation. 

Another case of great interest is that of an Aharonov- 
Bohm flux line in an integrable billiard p4fl. We show 
in Fig. 6 current streamlines of a state for the case of 
a unit square, a = b = 1, with a flux line at its cen- 
ter. The potential is this case is V(x) = — C, \x\ < -|, 
V(x) — +C, \x — 1| < |j repeated with period 2. Here 

C = (2n(/>/(j)o) k/VE. In Fig. 6 we have taken 0/0o = 
0.1, the flux in units of the flux quantum, to be small 
in order to minimize diffraction effects. The state is not 
much localized spatially, but the current is very regular. 
For symmetry reasons, the structure is particularly sim- 
ple along a coordinate axis, as we display in the figure. 
Because of diffraction, our theory of the flux line case 
has relatively large corrections. Nevertheless, our simple 
theory captures the main features of many states. We 
have related results J25| published elsewhere. 

In this case of a perturbed rectangular billiard, we 
have been able to find a bigger problem which admits 
a decomposition into fast and slow variables and thus 
a Born-Oppcnhcimer Ansatz. The original problem is 
solved by a superposition of results of the auxiliary prob- 
lem. 



VII. FAST AND SLOW ASYMPTOTICS 

We now give an example in which a separation of vari- 
ables into fast and slow does not hold over the whole 
domain, but does work over a limited domain, that is, 
however, sufficient to solve the problem. 

Consider a weakly distorted unit circle billiard. This 
is defined in polar coordinates by the boundary r = 1 + 
eAR(8), where e is small, AR is of order unity, and Ai? 
does not vary too rapidly. A number of papers have 
appeared on this topic recently M . 



?5i 






CM 
. 1 : 




H 


.3 - 






4 A 


"""vJXJ^^^ ^ 


: \ 


v. x_ — — ^_ 


c ~- 






1 1 1 1 1 1 1 1 1 1 1 1 1 


I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I. 



$ vanish at r = 1 + eAR(6) requires 



-.5 -.4 -.3 -.2 
x 



1 



k + keAR{9) 



v{ef 



2k 



a(9) 



(n--)Tr. (17) 



It is necessary to determine two functions, v and 
a, from Eq. jl7|) . We expect v to depend on 
\ [AR(0) + AR{9 ± tt)] = AR(6), and not on 5R(9) = 
4 [Ai?(0) — AR(9 ± 7r)] . Classically, a low angular mo- 
mentum state sees equally both sides of the circle, if at 
9, then also at 9 ± ir. To achieve this, we take a (9) = 
—keSR(9) + otQ, where cvo is to be determined. This im- 
plies that 



<e? 



E„ 



V{9) 



(18) 



where V{9) = 2k 2 eAR(9). This acts as a potential in the 
slow equation, and it will be important if k^/e >> 1. 

Again there is a Schrodinger equation with a periodic 
potential V{9) = V(9 + n), and thus solutions labelled 
by the band index m and 'crystal momentum' (3. Since 
physically ip(9 + 2n) = ip(9),(3 = or tt. The constant 
ao is determined from the special case AR = 0, which 
implies E m = m 2 , with m integer, (3 = tt (mmod2) . For 
= 0, the solution <fr(r\9) — J m (kr) implies ao — 
it/4. 

In Fig. 7 we show a state for a stadium billiard as in 
Sec. Ill, with a short straight side of length 2a — 2e — 
0.02 and end radius R=l. Then AR = e \cos9\ + 0(e 2 ). 
Some other states are shown in Rcf . S . 
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FIG. 6. Upper figure: Numerical persistent current stream- 
lines of a state induced by an Aharonov-Bohm flux line con- 
taining 0.1 flux quanta, located at the center of a square bil- 
liard. A quarter of the billiard is shown. Lower figure: Cur- 
rent in the y-direction on the y axis. The theoretical current 
is plotted. The numerical current is similar but has minima 
at the depths shown by the horizontal lines. The quantum 
numbers are n — 82, m — 0. 

We confine attention to high energy states whose clas- 
sical counterparts pass close to the center of the billiard. 
Such states have low angular momentum, that is, dimen- 
sionless angular momenta v << k. This suggests that we 
consider the radial variable r to be 'fast', and the angular 
variable 9 to be 'slow' and thus the Ansatz <I> = J„(0)(A;r). 
This, however, does not work. The reason is that near 
the origin, r is not fast compared with 9. 

However, we need to know $ only near the boundary 
in order to impose the Dirichlct conditions. Without a 
condition at the origin, using the asymptotic expansion 
for Bessel's functions suggests the Ansatz 



$(r|i 



1 



kr 



kr - 



v{9f 



2kr 



a{9) 



(16) 



where v{9) is the order of the Bessel functions and a(9) is 
a phase which mixes the asymptotic Bessel and Neumann 
functions. The slow equation is ip" + v 2 ip = 0. To make 




FIG. 7. Numerically obtained contour plot of a state of a 
Bunimovitch stadium billiard with end cap radius unity. whose 
straight sides are short, of length 2a, where a = 0.01. The par- 
allel lines near the left side of the boundary have separation 
2a. This state has m = 3, i.e. two transverse nodes, and 
n = 25. 

This application of the BOA differs from others in that 
we do not have a separation into slow and fast variables 
over the whole system, but we do have such a separation, 



asymptotically, over the part of the system that matters 
most. 



VIII. WHISPERING GALLERY MODES 

Quantization of higher order periodic orbit resonances 
in nearly circular billiards cannot be done easily by a 
Born-Oppenheimer method. However, the whispering 
gallery limit can be handled with a modification of the 
BOA, and is not restricted to the nearly circular case. 
Rather, we assume only a two dimensional smooth con- 
vex billiard. 

The results have long been known, but the methods 
we introduce are easier and more intuitive than earlier 
techniques. These modes correspond to classical motion 
which stays close to the billiard boundary while rapidly 
moving along the boundary. The effect was discussed by 
Lord Rayleigh, first for rays in his Theory of Sound, and 
later in terms of waves jll| . Lord Rayleigh did not con- 
sider the eigenmodes, however, which were first quantized 
by Keller |LJ] as an example of what came to be known 
as EBK quantization. Keller's work was based on the as- 
sumption of the existence of caustics in the corresponding 
classical motion. This was proved later by Lazutkin p6| ] 
who obtained an adiabatic invariant. 

We first obtain the result by SSTO. The surface of 
section is the billiard boundary. Let s be the distance 
along the boundary. The ray goes from s' to s, with 
s* = s — s' small. 

Let the radius of curvature at s = ^(s + s') be R(s). 
The curvature varies slowly in the sense dR/ds never gets 
too large and R is always positive and finite. Then, we 
use Eq. (0) above and approximate 



L(s,s') ~ 2i?(s)sin- 



s — s 



R(s) 
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We now look for a solution of Eq. (Q). We take as 
Ansatz 



ijj(s) = exp{ik [s — F(s)]} 



(20) 



where dF/ds = f(s) << 1. Thus, apart from the explicit 
exp(ifcs), tp(s) varies relatively slowly. 

We do the integral Tip by stationary phase. [The BOA 
below improves on this.] The stationary phase point is 
s' = s — s* where (s*) 2 /8R 2 = /. It can be checked later 
that the s' dependence of R can be neglected. Doing the 
integral gives 



Tip = i exp 
A solution requires 

f = fm{s)=Fi n 



i A -kR{2ff' 2 



3tt 



iP(s). 
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(22) 



This is the first quantization condition. The condition 
that / be small requires m « kR(s). 

The other quantization condition is a result of the re- 
quirement that ip is single valued on the boundary, that 
is, ip(0) = ip(£) where C is the circumference of the bil- 
liard. This condition can be expressed 



k [C - F m (C) + F m (0)] = 2irn 



(23) 



Thus the quantization depends on L R(s)~ 2 ^ 3 ds, effec- 
tively the mean 2/3 power of the curvature. 

We now turn to the BOA. The billiard is locally a 
circle of radius R(s). If R is constant, the solution is \P = 
J y {kr)e lu6 where J v is a Bessel function. We assume v is 
large. [This is close to the treatment of Lord Rayleigh, 
who, however, takes for granted that the index v is an 
integer. Most of the rest of his discussion concerns the 
asymptotic properties of Bessel functions of large index, 
a subject still under development at the time. The etalon 
method also considers this Bessel function, but does not 
distinguish fast and slow variables.] 

We assume that for a more general convex billiard W 
has this form locally. Namely we make the Ansatz 



#(r) = $(p\s)ip(s) = J„ {s) (kr s )ip(s) 



(24) 



where r s is a radial coordinate from the local center of 
curvature. Let the variable p > be given by r s — 
R(s) — p. With the exception of the factor exp(ifcs) all s 
dependence is slow compared with the 'fast' variable p. 
The index v can vary smoothly and there is no reason 
to make it an integer. We could already guess ip(s) = 
exp [i f s dsv(s)/R(s)] where we have replaced the local 
angle variable d9 s by ds/R. Parametrize v — kR(l — 
/) where / is small (and turns out to be the same / 
introduced earlier). Using a standard asymptotic formula 
|E7| we have 



(19) Jv(s)(k[R(s)- p]) = J v 
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We can replace v 1 ^ by (kR) 1 / 3 when multiplied with the 
small quantities / or p. 

The first quantization condition is that v(s) must be 
chosen to solve J v c s ) [kR(s)] — 0. A convenient analytic 
approximation is obtained from Eq. tea ) , again in terms 
of the zeroes z m of the Airy function, 



f = z m /2HkR)t. 



(26) 



This improves on Keller's treatment or the T operator 
above which is equivalent to approximating Ai(— x) oc 
sin(|:z;2 + 2L), valid for large x. 



To find ip(s) we substitute the Ansatz of Eq. (24) into 
the Helmholtz equation, neglecting derivatives of J v with 
respect to s. The angular derivative term (l/r s ) 2 d 2 /d9 2 
can be replaced by d 2 /ds 2 so we see that ip"(s) + 
(v 2 /R 2 )ip = 0. This has the approximate solution given 



in Eq. ( j20| ) above. From this, the requirement that the 
wavefunction be single valued on the boundary gives the 
second quantization condition. Thus the solution is 



V(s,p) ksC(s) Ai -z m + 



kp 



[kR(s)Y- 



i>(a) 



(27) 



The slowly varying prefactor C(s) is a normalization of 
the Airy function determined by current conservation. 

The validity of the BOA requires that the argument 
of the Airy function vary more rapidly with p than with 
s [in the region that the Airy function is large]. This 
condition can be written \R'\ << (kR) ' / z m . The ex- 
istence of a classical constant replacing fc 2 / 3 / ' z m in this 
strong inequality is the condition for the existence of a 
caustic. The caustic is at the zero argument of the Airy 
function, i.e. at p = z m (fci?) 1 / 3 /k. Thus if R is finite and 
smooth, there is always a sufficiently large energy such 
that states remaining close to the boundary exist. 

If 1/ R vanishes at some point, it is known that caustics 
do not exist, and the BOA must fail near that point. 
However, numerical evidence (at relatively low energies, 
of course) suggests that whispering gallery states exist 
even when the curvature vanishes. No detailed theory 
has been advanced for this case, to our knowledge. We 
have made some progress on this problem, which will be 
published elsewhere. 



IX. CONCLUSIONS 

With slight extensions, the Born-Oppenheimer ap- 
proximation can, rather quickly and easily, give the lead- 
ing order results for a number of interesting states. These 
include the important and well-known cases of the whis- 
pering gallery modes, and the bouncing ball states. It 
also includes some states only recently uncovered by the 
SSTO method. It is relatively easy to apply when, at 
least asymptotically, the variables can be separated into 
fast and slow. This can happen because of a small param- 
eter, or just because the particular state has that special 
property. 

The PEM and etalon methods are similar to the BOA 
in the sense that they deal directly with the partial dif- 
ferential equations. The detailed procedures and moti- 
vations for each step is rather different from the BOA. 
The leading results are the same, however. An advantage 
of the BOA is that it is simpler and taught in standard 
courses on quantum mechanics while the other methods 
are less well known. 

These methods are not a priori semiclassical. How- 
ever, if a variable is fast because it has more energy than 
the slow variable, semiclassical methods can be used. The 
Keller-Rubinow method is classic, but a little cumber- 
some in practice. The SSTO method seems to be in a 
certain sense more general than the Born-Oppenheimer. 
Compared with Keller-Rubinow, it has the advantage of 
the simplifications coming from use of the surface of sec- 
tion. Compared with BOA, it first finds the slow wave- 



function, and from that the fast one. Often the slow 
wavefunction is more interesting physically. 

Another well known semiclassical method sometimes 
compared with the SSTO method is the Birkhoff- 
Gustavson normal form |2S|| . However, Birkhoff- 
Gustavson is not adapted to finding special wavefunc- 
tions. It rather gives a large number of wavefunctions 
and energy levels at once, and is similar to a numerical 
method in that respect. 

Formal corrections to the PEM and etalon methods 
have been written down p7| , which differ from the usual 
methods used with the BOA. These corrections are not 
much used in practice. Rather, one resorts to a numerical 
method. Moreover, quite large corrections to the wave- 
functions are sometimes found in the numerics. This is 
because the theory singles out some relatively small sub- 
set of special states which are approximately decoupled 
from all the rest. Nothing prevents some unrelated state 
from having nearly the same energy as the special state. 
The residual coupling then mixes the two appreciably. 
The energies are given quite well, although not neces- 
sarily as accurately as the mean level spacing of all the 
levels. 

Although it is sometimes possible to find approxima- 
tions to both states and to estimate the mixing param- 
eters, it is hard to do that systematically. It is perhaps 
more advisable to think of the wavefunctions obtained 
as being those of quasimodes p9j| , in other words, as lin- 
ear superpositions of a few nearly degenerate true eigen- 
states. In many situations, of course, a quasimode pro- 
vides a correct and more physical description of a phe- 
nomenon than the true eigenmodes. 
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